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Introduction 

This document describes the Discrete Fourier Transform (DFT), that is, a Fourier Transform as applied 
to a discrete complex valued series. The mathematics will be given and source code (written in the C 
programming language) is provided in the appendices. 

Theory 

Continuous 

For a continuous function of one variable f(t) 5 the Fourier Transform F(f) will be defined as: 



00 



F(f)= \ 



f(t) e"j 2Tm dt 



and the inverse transform as 



00 




F(f) ei 2flft df 



-oo 



where j is the square root of -1 and e denotes the natural exponent 



eJ 0 = cos(tf) + j sln(0). 



Discrete 



Consider a complex series x(k) with N samples of the form 



XO, X1, *2, X 3 ■ ■ - X k . . . XfM 



where x is a complex number 
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Xj - X r eal + J x imaq 

Further, assume that that the series outside the range 0, N-1 is extended N-periodic, that is, x k = x k+N for 

all k. The FT of this series will be denoted X(k), it will also have N samples. The forward transform will 
be defined as 

N-1 

X(n) = 77 2 x(k) e ~i k2T7n/N for n=O..N-1 
IN k=0 

The inverse transform will be defined as 

N-1 

x(n)= 2 X(k) e j k27Tn/N forn=O..N-1 
k=0 

Of course although the functions here are described as complex series, real valued series can be 
represented by setting the imaginary part to 0. In general, the transform into the frequency domain will 
be a complex valued function, that is, with magnitude and phase. 

magnitude = ||X(n)|| = ( X re al*Xr ea l + Ximaq*Ximaq ) 0 5 

phase = ten- 



The following diagrams show the relationship between the series index and the frequency domain 
. sample index. Note the functions here are only diagramatic, in general they are both complex valued 
series. 




Notes 



The first sample X(0) of the transformed series is the DC component, more commonly known as 
the average of the input series. 
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• The DFT of a real series, ie: imaginary part of x(k) = 0, results in a symmetric series about the 
Nyquist frequency. The negative frequency samples are also the inverse of the positive frequency 
samples. 

• The highest positive (or negative) frequency sample is called the Nyquist frequency. This is the 
highest frequency component that should exist in the input series for the DFT to yield 
"uncorrupted" results. More specifically if there are no frequencies above Nyquist the original 
signal can be exactly reconstructed from the samples. 

• The relationship between the harmonics returns by the DFT and the periodic component in the 
time domain is illustrated below. 



DC 
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DFT and FFT algorithm. 
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While the DFT transform above can be applied to any complex valued series, in practice for large series 
it can take considerable time to compute, the time taken being proportional to the square of the number 
on points in the series. A much faster algorithm has been developed by Cooley and Tukey around 1965 
called the FFT (Fast Fourier Transform). The only requirement of the the most popular implementation 
of this algorithm (Radix-2 Cooley -Tukey) is that the number of points in the series be a power of 2. The 
computing time for the radix-2 FFT is proportional to 



So for example a transform on 1024 points using the DFT takes 10 times longer than using the FFT, a 
significant speed increase. Note that in reality comparing speeds of various FFT routines is problematic, 
many of the reported timings have more to do with specific coding methods and their relationship to the 
hardware and operating system. 

Sample transform pairs and relationships 

• The Fourier transform is linear, that is 



N log2(N) 



af(t) + bg(t) 



->aF(f) + bG(f) 



ax k + by k 



->aX k + bY k 



• Scaling relationship 



f(t/a) 



-> a F(a f) 



f(at)--> F(f/a)/a 



• Shifting 



f(t + a)~- > F(f) e-j 2 P { 



• Modulation 



f(t)ei 2 P iat -> 



F(t - a) 



• Duality 



X k ->(l/N)x N . k 



Applying the DFT twice results in a scaled, time reversed version of the original series. 



• The transform of a constant function is a DC value only. 



• The transform of a delta function is a constant 
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4 

in L 



The transform of an infinite train of delta functions spaced by T is an infinite train of delta 
functions spaced by 1/T. 



1/T 



The transform of a cos function is a positive delta at the appropriate positive and negative 
frequency. 




4 

1/2! 



The transform of a sin function is a negative complex delta function at the appropriate positive 
frequency and a negative complex delta at the appropriate negative frequency. 

1/2 i 




-1/2 



The transform of a square pulse is a sine function 



w 




1/W 



More precisely, if f(t) = 1 for |t| < 0.5, and f(t) = 0 otherwise then F(f) = sin(pi f) / (pi f) 
• Convolution 

f(t)xg(t)->F(f)G(f) 
F(f) x G(f) -> f(t) g(t) 
x k xy k ->NX k Y k 

x k y k ->(i/N)X k xY k 
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• Multiplication in one domain is equivalent to convolution in the other domain and visa versa. For 
example the transform of a truncated sin function are two delta functions convolved with a sine 
function, a truncated sin function is a sin function multiplied by a square pulse. 

o The transform of a triangular pulse is a sine 2 function. This can be derived from first principles 
but is more easily derived by describing the triangular pulse as the convolution of two square 
pulses and using the convolution-multiplication relationship of the Fourier Transform. 

Sampling theoreno 

The sampling theorem (often called "Shannons Sampling Theorem") states that a continuous signal must 
be discretely sampled at least twice the frequency of the highest frequency in the signal. 

More precisely, a continuous function f(t) is completely defined by samples every 1/fs (fs is the sample 
frequency) if the frequency spectrum F(f) is zero for f > fs/2. fs/2 is called the Nyquist frequency and 
places the limit on the minimum sampling frequency when digitising a continuous sugnal. 

If x(k) are the samples of f(t) every 1/fs then f(t) can be exactly reconstructed from these samples, if the 
sampling theorem has been satisfied, by 

k=oo 

f<0= I x(k) slncd fs - k) 
k=-oo 



where 



sinc(x) 



sin(TJx) 

TtX 



Normally the signal to be digitised would be appropriately filtered before sampling to remove higher 
frequency components. If the sampling frequency is not high enough the high frequency components 
will wrap around and appear in other locations in the discrete spectrum, thus corrupting it. 

The key features and consequences of sampling a continuous signal can be shown graphically as 
follows. 

Consider a continuous signal in the time and frequency domain. 





Sample this signal with a sampling frequency f g , time between samples is l/f g . This is equivalent to 
convolving in the frequency domain by delta function train with a spacing of f g . 
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FT 



kflL 





fs 



0 fs fs 
2 



If the sampling frequency is too low the frequency spectrum overlaps, and become corrupted. 




/ w\ 



0 fs fs 
2 



Another way to look at this is to consider a sine function sampled twice per period (Nyquist rate). There 
are other sinusoid functions of higher frequencies that would give exactly the same samples and thus 
can't be distinguished from the frequency of the original sinusoid. 
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